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ABSTRACT 

In this paper the alternate forms of the relative attitude kinematics and relative dynamics equations are 
presented. These developments are different from the earlier developments that have been presented in other 
publications. The current forms of equations have the advantage of being simpler than earlier ones. These equations 
are applied in developing the necessary kinematics and dynamics for relative navigation in formation flying and 
virtual platforms. These equations also have application in the implementation of nonlinear full state feedback and 
nonlinear output feedback control for large attitude angle acquisition and tracking. This paper presents simulations 
from such a full state feedback control application. 

INTRODUCTION 

Navigation and control for spacecraft flying in formation requires the relative attitude information. Since the 
attitude rotation matrices and their kinematics relationships are generally defined in the inertial reference frame, 
Ref.l was the first attempt at completely defining the kinematics and relative attitude dynamics when the attitude 
matrix is defined in a non-inertial frame. While the concept of relative attitude had earlier been addressed in Refs.2- 
4, Ref.l was the first publication to provide the complete development that is necessary for formation flying control. 
This paper continues these developments, providing the relative attitude kinematics and dynamics that are 
alternative to the forms developed in Ref.l. 

The attitude of a rigid-body with respect to the inertial frame is determined by a rotation transformation matrix 
from the inertial frame to body frame. This rotation matrix is referred to as the attitude matrix. In practical design, as 
noted in Ref.l, the attitude matrix is parameterized to be 4-dimension parameters such as axis/angle variables and 
quaternion, and 3-dimension parameters such as Euler angles, Rodrigues (Gibbs vector) and modified Rodrigues,. 
Similarly, rotation matrix of the body-frame with respect to a non-inertial frame can be defined as the relative 
attitude matrix. 

In developing the kinematics and dynamics of the relative attitude, Ref.l addressed the following: 

i. Kinematics equations for relative attitude matrix 

ii. Kinematics equations for relative attitude parameters 

iii. Relative attitude dynamics equations 

This paper proceeds along the same line, providing alternate solutions to the above items. The primary 
difference in this approach results from a different definition of the relative angular velocity. This alternate form 
provides simpler relative kinematics and dynamics equations than those in Ref.l . 

SYMBOLS 

[dp 1 , dp2, dp3] the three components of relative modified Rodrigues parameters 
[dwl, dw2, dw3] the three components of the relative angular velocity ( rad/sec) 
g Rodrigues parameter 

p modified Rodrigues parameter 
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q quaternion parameter 

R bd relative attitude matrix (transformation matrix) of reference frame with respect to .y D . 

co angular rate 

d(.)/dt B vector time derivative in 

d(.)/dt D vector derivative in :rF n 

vectrice of the inertial reference system 
._?" x vectrice of the local horizontal reference system for target satellite 

vectrice of the local horizontal reference system of the chasing spacecraft 
•?" B vectrice of the body-fixed reference system of the chasing spacecraft 

&~ D vectrice of the body-fixed reference system of the target attitude 

<8> nonlinear operator 

REPRESENTATION OF RELATIVE ATTITUDE 


Relative Attitude Matrix 

It is assumed that any reference frame can be denoted by a vectrix 7 . The following series of reference frames 
can then be defined. 
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The attitude of a rigid spacecraft is the orientation of the reference frame with respect to another frame. The 
most convenient reference frame is a dextral, orthogonal triad which is fixed with the rigid body of spacecraft. The 
other reference frame can be an inertial reference frame, or it can be a moveable reference frame which is fixed to 
another body. The attitude with respect to the inertial reference frame is the absolute attitude. The attitude with 
respect to a movable rotting reference frame is the relative attitude, which is respect to a movable rotating reference 
frame is named the relative attitude. For example in order to study the orientation relationship between two rotating 
reference frames, the relative attitude can be defined as a transformation matrix between two rotating reference 
frames. 
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Parametric Representation of the Relative Attitude Matrix 

As with parametric representation of absolute attitude matrix (Refs. 6-8), the relative attitude matrix can also be 
represented by the attitude parameters such as Axis/angle, Quaternion, Rodrigues parameters (Gibbs vector), 
modified Rodrigues parameters and Euler angles. The relationship between relative attitude matrix and absolute 
attitude matrix is 


^BD^Pbd) Rb£Pi) R DliPd ) 


(3) 


It is obvious that the relative attitude parameters g bd and absolute attitude parameters g b , g d have a nonlinear 
relationship. This can be represented by the following unified notation: 

(4) 


Ptd = P b ®Pd 
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where the rule of the nonlinear operator ® is determined by composition rotation rule of the attitude parameters [6], 
This operation will, therefore, be different for different attitude parameters. 


For the modified Rodrigues parameters: 

Pu . Pt ® Pi ' - pjpJpr'-yp^-plpj-^p: k 

1 + pJPdPb r Pb +2 PJPb 


For Rodigues parameters (Gibbs vector): 



g b -g d + \sl]g d 
1 + SbSd 


For Quaternion parameters: 
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RELATIVE ATTITUDE KINEMATICS FOR ATTITUDE MATRIX 

Equation (2) can be written as 



Taking derivative with respect to the inertial reference frame for both sides, 




x b R bd) 




+ co x 
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Based on the differential rules for the vectrices (Ref.7), Eq.(10) can be developed as 
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( 11 ) 


dR„ 


<oX^ b ^ T bR bd =X-~ + 
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In the body reference frame ,^ B , the kinematics equation for relative attitude matrix R BD will be 


dR B 


dt„ 


]i? SD C0 A ) ]R BD 


( 12 ) 


If relative angular velocity is defined as 


ffl. , = m - £0 , 

bd b d 


X( a b- R BD a d) ®L <^ b - R BD^d) 


( 13 ) 


the attitude kinematics equation for the relative attitude matrix (12) can be written the compact form in the body 


reference frame 
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Considering 


^BD^DB I 


( 15 ) 


the relative kinematics equation for R DB in Body reference frame .j^ B becomes 

d R DB _ r> r , b \ x-» 


dt 


( 16 ) 


In the similar way, the kinematics equation for R DB in reference frame J# D can be obtained: 
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where 
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(19) 
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these kinematics equations for relative attitude above can be represented in form of a theorem of relative attitude 
kinematics: 

Theorem of Relative Attitude Kinematics 

The angular velocities of with respect to the inertial reference frame are are co b ,co d , and 


= ^ B o h 


= K®d 


( 20 ) 


The relative attitude between the reference frames .y , , , ,3~ D are: 

B ~ 1 ^~D ~ ^Df^l ^B = ^BD’-^D 


( 21 ) 


The kinematics equations for R BD and R DB are: 

Kd = -[(coL) x ]^=[(4) X ]^ D =^n[(4) x ] = -R BD [(<dT\ (22) 

r db = -K DB \-(K b r\ -R DB [(Kn =«</] - -^J(4n (23) 

where 


(®i ® db R B D®db ^BD^bd (24) 

RELATIVE ATTITUDE KINEMATICS FOR ATTITUDE PARAMETERS 

The First Kind of Relative Attitude Kinematics Equations 

Based on the definitions of the relative attitude and relative angular velocity given above the relative attitude 
kinematics equations for various attitude representation have been developed in Ref.l . These equations are referred 
to as the first kind of relative attitude kinematics equations, and are summarized as follows: 

Quaternion 

Relative kinematics equation 


Vbd = |([9mWmA°L 




(25) 
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Relative velocity 


®bd -(QbdAbd P-bdflbd Qb^^bd) 


( 26 ) 


Attitude transformation matrix 

R BD = ( C Ibd4- < lb<flbd)I + 2 e Jbcflbd~2qbd4[ c Ibd X ] 


( 27 ) 


Rodrigues Parameters 

Relative attitude kinematics equations 

Sbd = |(kw] + S b Sbd + I)*bd 


Relative velocity 



— -r-v-isldMbd 

(! + SbcSbd) 


( 28 ) 


( 29 ) 


Attitude transformation matrix 

r bd = - — x ^:^^Jg bd )i^g bSbd T -2[g b /]) 

Modified Rodrigues Parameter 
Relative attitude kinematics equation 

Pbd = ^K 1 -PbcPbd) I+2 \Pl2 +2 PbcPbd} b hd 

Relative velocity 

®w = 7 -K 1 -PkdPj+lPbcPbd-ZiPbd]! P bd 

(1 +pLp,,) 2 


Attitude transformation matrix 


R BD = 7 -[( l - 6 PbdPbd + (Pb T t Pbd) 2 ) 1 + ZPbcPL-^ -PbdPbd)\Pbd\] 

0 + PbcPbd) 
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The Second Kind of Relative Attitude Kinematics Equations 

The relative kinematics equations are dictated by the definition of the relative angular velocity. Using the 
following (alternate) definition of angular velocity 

Aco=co A -co rf 


an alternate set of relative attitude kinematics are developed. These may be defined as the alternate (second) form of 
relative attitude kinematics equations. It should be noted that in this case, the subscript of relative attitude are 
replaced by the use of A. 

Quaternion 

Using the new notation, the first kind of relative kinematics equation in Quaternion can be written as 

*** =m4 

dt bd 


(35) 


where 
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the relative attitude kinematics Eq.(35) can be written as 

^ =MAco +M(I-R BD )o) d 


(38) 


where 

R bd = {Aqj-Aq r Aq)I+2AqAq T -2Aq 4 [Aqx] (39) 


Because 
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Considering the following equality relationships for any 3 dimensional vector p: 


( 40 ) 
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\P*]\P*] + P T PI-PP T = 0 3 X3 , 


(41) 


T T T T 

p ppp -pp pp 


PP T \P *] =0 3 x3’ \P r=0 3x3 


Eq.(40) can be simplified as follows: 
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Finally Eq.(38) can be simplified to obtain the second kind of relative attitude kinematics equations 
in Quaternion : 
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Relative angular velocity: 


Aco = 2(Aq 4 Aq -Aq.Aq -AqxAq) +{R BD (Aq) -7)co d (45) 

Attitude transformation matrix: 

R BD (Aq) = (AqJ-Aq r Aq)I+2AqAq T -2Aq 4 [Aqx] (46) 


Rodrigues Parameters 

Using new notation, the relative attitude kinematics equations can be written as 

~ -MAw +M(I-R BD )(o d (47 ) 


where 
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Applying equality relationship Eq.(41-42) simplifies the equation as 

W~R bd > ~([4? X ] + A^A g r + /)(/- l- ((1 -Ag T Ag)I + 2AgAg T -2[Ag'})) 

2 (l+A^Ag) 


= [Ag*] 


( 49 ) 


Substituting Eq.(49) into Eq.(47), the second kind of relative attitude kinematics equation in Rodrigues 
parameters are obtained: 


~ = -|([Ag x ] + AgAg T + 7)A(o+[Ag >< ]co rf 


( 50 ) 


Relative velocity: 


A® = [Ag x ])Ag+(R BD (Ag)-l)ca d 

(1+Ag'Ag) 


( 51 ) 


Attitude transformation matrix: 

R BD (Ag) = l— {(l-Ag T Ag)I + 2AgAg T -2[Agx]) 

(l+A^Ag) 


( 52 ) 


Modified Rodrigues Parameter 

In a similar way, the second kind of relative attitude kinematics equations in modified Rodrigues parameters 
can be written as 


^ = — [(1 -A p r Ap)I+2[Ap x ] +2A pAp r ]A©+[Ap *]co rf 


( 53 ) 


Relative velocity: 


Aca = — -[(1 ~Ap T Ap)+2ApAp T -2[Ap *]] Ap+[R BD (Ap)-I\(d d 

(1 +Ap t A pf 


( 54 ) 
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Attitude transformation matrix: 


R bd( A p) = ~ rt(l -6A p T Ap+(Ap T Apf)I + SApAp r -4(l -A p T Ap)[Ap *]] ( 55 ) 

(1 +Ap Ap) v ' 

RELATIVE ATTITUDE DYNAMICS EQUATIONS 

The First Kind of Relative Attitude Dynamics Equations 

It is assumed that J b (t), co b (t), h b (t) and L,, are the inertial dyadics, angular velocity, angular momentum and the 
applied angular moment of the chase satellite. Similarly, J,j(t) , co^t) , h^t) and Lj are the corresponding items for 
the target satellite. The relative angular velocity of the chase satellite body reference frame with respect to the 
target satellite reference frame is defined as follows: 


m w =£ v^ = 




(56) 


The relative attitude dynamics equation in pursuer satellite body reference frame developed in Ref.l is 

'V* ) M + [( < 0 m) \-h (l ^hd + i( R BD Cj) d) Vb 0} bd + [(®m) Vb^BD^d^b^^bd) 

~ R BD^d~^BD^ A ^cPcl + \- (Si JAJ/0^) 


(57) 


where 


+ KVA = L d AJ d = R DB J b R BD ~J d 


(58) 


The Second Kind of Relative Attitude Dynamics Equations 

Using the alternate definition of relative angular velocity: 

Ao)=ffl,-ra rf ( 59 ) 


The relationship between co bd and Aco (the two definitions of relative angular velocity) is 

®L =Aro+(/-JVco rf (60) 

Considering Eq.(59) and (60), the left hand terms of Eq.(57) can be developed as: 

A^fcCTiAfb "^(/-Rg^CbUJ^AO) 1 R BD^d +R bP^c^BD^d (61) 
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[(©L)1 J i, ra M = [ A ®1 J A A “ + [®^4 A ®“[(^BD®rf)1 J 6 A ® + [ A ®l4 C0 d't A ® X ] J / J SD C0 rf 

T [®d^®rf _ [®d^^a/)®d _ [(-^g£i®^) ]^4®rf + [(^B£)®rf) Vb^BD^d 


( 62 ) 


[(®i<f) Vb^BD^d t A ® ^b^BD (£> d + ^aV h^BD^d K^Bi)®rf) Vb^BD^d 


(63) 


[(^BD®d) Vb a bd tC^B£>®rf) Vb A ® + [(^B£)®d) Vb^d K^Bn®rf) Vb^BD^d 


(64) 


^6[(®w) 1 ^B£>®d /f,[ A ® -I ^BD®B + ^b t ®rf] -^fiD® rf At(^BD®rf) ]^BD®rf 


(65) 


Substituting Eqs. (61-65) into Eq.(57) and re-arranging 

4 A ® + [Acd x ]./ 6 A (0 + [Aco*] J b (O d = J b (R BD -I)(z d - [co> a Acd - [co^] J b a d 
+ K^flO®^ ^b^BD (l} d + ^b “ ^BD^d~^BD^^/^d + [®rfl A 4 /®tf) 


( 66 ) 


where 


J >rf + 


(67) 


These are the second kind of relative attitude dynamics equation. The two equations may be summarized in terms of 
a relative attitude dynamics theorem: 

Theorem of Relative Attitude Dynamics 

If the relative angular velocity is defined as 


®M=® A -^n® rf 


( 68 ) 


then the first kind of relative attitude dynamics equation in the chase satellite body reference frame is 

4®L + [(®L) x ]A®L + [(®L) x ]7An®, 

= 7 6 [(C0 m ) ] R-BD^d ~ [(*B 0 ®B) Vb a bd + ^b~^BD^d~^BD^^J S> d + ^ (i) d\^cf 0 e) 


(69) 
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If the relative angular velocity is defined as 


Aa=a b -(S} d 


( 70 ) 


then the second kind of relative attitude dynamics equation in the frame is 

J 6 Affl+[Aco x ]J i Affl+[Aoj x ]y 4 co rf =J 6 (i? 5 Z) - J 06 d -[toXAra-[co^J 4 (o rf 

BD a d + ^b ^BD^d R BD (UJi> d +[<Sb d \A.J t p^ 


( 71 ) 


where 


^ d ^DB^b^BD J d 


( 72 ) 


J Pd + = L d 


( 73 ) 


APPLICATION TO LARGE ATTITUDE ANGLE ACQUISITION AND TRACKING CONTROL 

An application of relative attitude kinematics and dynamics equations developed in this paper is the large angle 
acquisition and tracking control. This section provides the simulations of such an application. The advantage in 
using relative attitude is that the tracking control problem is converted into a regulator problem, simplifying the 
control system design. 

It is assumed that the second kind of relative kinematics and dynamics equations in modified Rodrigues 
parameters are used for the acquisition and tracking control. The equations are as follows: 

Relative kinematics equations: 

= MA®+[Ap*](i) d (74) 


where 


M = — [(1 -Ap r Ap)I+2[Ap*]+2ApAp T ] 
4 


( 75 ) 


The relative dynamics equations: 


JA&> + Acox M d + Acox JAco =L 


( 76 ) 



where 


1 =J b( R BD - f ^d ~ K'14 A “ - [®X®d 

[(R BD a> d ) ].J h R BD (a d ^L b -R BD L d ~R BD (AJ i p d +[(ij d \kJ < P'^ 


( 77 ) 


The Lyapunov function V can be selected to be 


V=K p Ap T Ap+—A(o T JAa> 


( 78 ) 


where Kp is a positive constant. The time-derivative of Lyapunov function V is 

V=2K p Ap T Ap+Aa T JAa =Aa T (2K p M T Ap+L) 


( 79 ) 


If 


L = -2K p M T A p -KjJ Aco 


( 80 ) 


then 


V^-KjA co r JAcos 0 


( 81 ) 


where K d is a positive constant. The closed loop system dynamics equation is 

JA(a + Ao>x Jw d +A®x JAco = ~K (l JA^-2KM J Ap 


( 82 ) 


Equation (81) implies that V(t)< V(0), and therefore, that A p and A to are bounded. In addition, from Eq.(81) 

V=-2K d AmJA(s> ( 83 ) 


so it can be seen that d 2 V/dt 2 is bounded. Hence dV/dt is uniformly continuous [9], Application of Barbalat’s 
lemma [9] then indicates that Aco - 0 as t- °°. Considering the closed loop equation (82), it can be obtained that 
Ap-0 as t- °°. Therefore, we have that (Ap, Aco) - (0,0) as t - . It means that the nonlinear control law given by 

Eq.(84) 


R b R BD^d Jb( R BD •0®rf + [ C °J l 4 Ac0 + [ C0 (/]7 A (O rf 
“1C Ap~KjJAcd 


( 84 ) 


is a global asymptotically stable control law for the system given by Eqs. (74-77). 
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SIMULATIONS 


The Lyapunov nonlinear attitude control law has been used for large attitude angle acquisition and tracking 
control for the EO-l/LandSat 7 formation. For this simulation the control is implemented using full state feedback 
only. The simulations using measurement output feedback will be presented in future papers. 
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Figure 1. Large Attitude Angle Acquisition and Tracking Control for EO-l/LandSat7 
Formation using Full State Feedback 


The initial conditions for this full state feedback simulations are as follows: 


Maneuver angle = 70 deg. 

Target satellite (Landsat 7) angular velocity = (00 co x ) T 
Pursuer satellite (EO-1) angular velocity = ( 0.1 co T 0.01oo T 2co x ) T 
co T =angular velocity of the 705 km circular orbit 

[dpi, dp2, dp3] in Fig 1 are the three components of the relative attitude represented in the relative modified 
Rodrigues parameters, and [dwl, dw2, dw3] are the three components of the relative angular velocity ( rad/sec). In 
this simulation, it is assumed that the attitude state vector is given and there no attitude measurement information 
available for feedback control. 


CONCLUSIONS 

In this paper the alternate relative attitude kinematics and dynamics equations are developed for the various 
attitude parametric representations. Compared with the first kind of relative attitude kinematics and dynamics 
equations (Ref.l), this has the advantage of being simpler. These developments will find ready application in the 
problems of relative attitude determination and control, and will be very useful for spacecraft formation control and 
relative navigation. As an example of such application, the Lyapunov nonlinear control law for large attitude angle 
acquisition and tracking has been developed and simulated for the EO-l/LandSat 7 formation. This simulation 
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implemented the foil state feedback control. 
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